Materials Design using Correlated Oxides: Optical Properties of Vanadium Dioxide 
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Materials with strong electronic Coulomb interactions play an increasing role in modern mate- 
rials applications. "Thermochromic" systems, which exhibit thermally induced changes in their 
optical response, provide a particularly interesting case. The optical switching associated with the 
metal-insulator transition of vanadium dioxide (VO2), for example, has been proposed for use in 
"intelligent" windows, which selectively filter radiative heat in hot weather conditions. In this work, 
we develop the theoretical tools for describing such a behavior. Using a novel scheme for the cal- 
culation of the optical conductivity of correlated materials, we obtain quantitative agreement with 
experiments for both phases of VO2. On the example of an optimized energy-saving window setup, 
we further demonstrate that theoretical materials design has now come into reach, even for the 
particularly challenging class of correlated electron systems. 
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The concerted behavior of electrons in materials with 
strong electronic Coulomb interactions causes an extreme 
sensitivity to external stimuli, such as temperature, pres- 
sure or external fields. Heating insulating SmNi03 be- 
yond 400 K or applying a pressure of just a few kbar to 
the Mott insulator (Vi-aCr^Oa (x=0.01), e.g., make 
the materials undergo transitions to metallic states [J]. 
This tuncability of fundamental properties is both a 
harbinger for diverse technological applications and a 
challenge for a theoretical description. Indeed, many 
materials used in modern applications fall in this class 
of strongly correlated materials where electron-electron 
interactions profoundly modify (if not invalidate) a pure 
band picture. State-of-the-art first principles methods, 
such as density functional theory, are then no longer suf- 
ficient to predict physical properties of these compounds. 

An increasing role is nowadays played by artificial 
structures, ranging from spintronic multi-layers 0, 01 to 
functional surfaces, where appropriate coatings e.g. pro- 
vide a self-cleaning mechanism |J|. However, the huge 
freedom in the design, which concerns not only the chem- 
ical composition, the doping and the growth conditions, 
but also the geometry of the setup (e.g. the layer thick- 
nesses), makes the search for devices with specific elec- 
tronic properties a tedious task. Ideally, the quest for 
promising materials and setups should therefore be sec- 
onded by theoretical predictions. Yet, the abundance 
of fruitful experimental work on tailoring e.g. materials 
with specific optical absorption profiles is met with only 
scarce contributions from theoretical investigations Q . 

Here, we demonstrate that even for the particularly 
challenging class of correlated materials a quantitative 
description, and thus predictions, are possible. We in- 
troduce a novel scheme for the calculation of optical 
properties, which can in principle be employed in con- 
junction with any electronic structure technique that 
uses localized basis functions. We shall demonstrate the 



power of the approach within the framework of the com- 
bined density functional - dynamical mean field method 
"LDA+DMFT" ^ , on the example of the optical con- 
ductivity of vanadium dioxide, VO2. We further inves- 
tigate VO2 based heterostructures and show that the- 
oretical materials design of correlated material-derived 
devices is coming into reach. 

Within dynamical mean field theory, emphasis is com- 
monly put on spectral properties, and the evaluation of 
observables other than spectral functions is a rather new 
advancement in the realistic context. Yet, as explained 
above, it is the response behavior of correlated systems 
that is promising for applications. In this vein, recent pi- 
oneering work 0, d, 0, E3 has successfully described 
qualitative features of the optical response of several cor- 
related materials. The calculation of absolute values in 
theoretical response functions - a prerequisite for quanti- 
tative theoretical materials design - has, however, turned 
out to be a formidable challenge. This is mainly due to 
the sensitivity of the obtained results on the accuracy of 
the transition matrix elements 
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In the framework of linear response theory, and when 
neglecting vertex-corrections [Hi ] , the optical conductiv- 
ity can be expressed as (for reviews see [H, [13]) 



Rea af3 (uj) 



2ne 2 h 
V 



E 



c\lu 



, f(w') - f(u/ + LO) 



x tr<^ A k (cj' + uj)v Ka A k (uj')v k ^ \ (1) 



Correlation effects enter the calculation via the spec- 
tral functions A k (u>), while the Fermi velocities v ktCl = 
— (kL'I'PQ.Ik-L), matrix elements of the momentum oper- 
ator V, which weight the different transitions, are deter- 
mined by the band-structure. Both, spectral functions 
and velocities are matrices in the basis of localized or- 
bitals, indexed by L = (n, 2,771,7), with the usual quan- 
tum numbers (71, /, m), while 7 denotes the atoms in the 
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FIG. 1: Optical conductivity of (a) metallic, (b) insulating VO2 for polarizations E. Theory (red) ([aab] = [0.85 0.85 0.53]), 
experimental data (i) single crystals [f3| (green), (ii) thin film [lj (solid blue), (iii) poly crystalline film [T3 (dashed blue). 



unit cell : |kL) is the Fourier transform of the Wannier 
function x RI (r) localized at atom 7 in the unit cell R. 
The Fermi functions f(w) select the range of occupied 
and empty energies, respectively, V is the unit-cell vol- 
ume, a, (3 denote cartesian coordinates, and Re a af3 is the 
response in a— direction for a light polarization E along 
/?. Approaches such as LDA+DMFT [|| require the use 
of localized, Wannicr-likc, basis sets, rendering the eval- 
uation of the full matrix element tedious. To this end, we 
developed a "generalized Peierls" approach that extends 
the formalism for lattice models (see the review 17[) to 
the realistic case of multi-atomic unit cells. We find that 
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Here, pl denotes the position of an atom within the unit- 
cell. The term in brackets, which is used in the actual cal- 
culations, we refer to as the "generalized Peierls" term : 
While the derivative term is the common Fermi velocity, 
the one proportional to the Hamiltonian originates from 
the generalization to realistic multi-atomic unit-cells. 
The correction term that recovers the full matrix element 
is denoted T (for its explicit form see [Hj]). The latter 
reduces to purely atomic transitions, (R, 7) = (R',7'), 
in the limit of strongly localized orbitals. In other words, 
the accuracy of the approach is controlled by the choice of 
basis functions. This generalized Peierls approach is thus 
expected to be a. good approximation for systems with lo- 
calized orbitals, such as the 3d or 4f orbitals in transition 
metal or lanthanidc/actinidc compounds. This expecta- 
tion is indeed true as can be verified numerically [13]. 
From the optical conductivity, one can calculate the spec- 
ular reflectivity and transmittancc. Therewith, also the 
specular color of the material becomes accessible [34 1. 

Our focus compound, VO2, has attracted a lot of at- 
tention recently [l8|, [l9(. It is among those materials 
in which correlation effects play a decisive role, to the 
extent that standard band-structure approaches fail to 



capture even most basic experimental facts : In metallic 
VO2, important incoherent spectral features witnessed 
by photoemission experiments are absent in band the- 
ory, and the metal-insulator transition at T c =340 K is 
missed altogether. For reviews see e.g. [l2|, HO] and 
references therein. Hence, any description of the ther- 
mochromic properties of this material must foot on an 
electronic structure approach that goes beyond band- 
theory and masters the many-body effects at work. In- 
deed, LDA+DMFT results for the spectral properties of 
VO2 agree well with experimental findings in both, the 
metallic U 2^, [23| and the insulating phase [23[ . Foot- 



ing on the electronic structure calculation of Ref. [23[ 
and our recent extension thereof 24 , 25[ , we first employ 
our scheme to access optical properties of bulk VO2 both 
above and below its metal-insulator transition. 

Fig.QJa) shows our results for the optical conductivity 
of metallic VO2 as a function of frequ ency and in com- 
parison with experimental data [111, [ll, [TBI], see a l so G3- 
As can be expected from the crystal structure j2(j], the 
optical response depends only weakly on the light po- 
larization. The Drude-like metallic response is caused by 
transitions between narrow vanadium 3d excitations near 
the Fermi level, and thus affects only the low infra-red 
regime - a crucial observation as seen in the following. 
The shoulder at 1.75 eV still stems from intra- vanadium 
3d contributions, while transitions involving oxygen 2p 
orbitals set in at 2 eV, henceforth constituting the major 
spectral weight up to the highest energies of the calcula- 
tion. Also for insulating VO2, see Fig. [TJb), the results 
are in good agreement with experiments. This time, a 
slight polarization dependence is seen in both, experi- 
ment and theory, owing to the change in crystal sym- 
metry. Indeed, vanadium atoms pair up in the insulator 
to form dimers along the c-axis, leading to the forma- 
tion of bonding/anti-bonding states for 3d orbitals [28j |. 
Optical transitions between these orbitals result, in the 
corresponding energy range (w=1.5-2.5 eV), in a higher 
amplitude of the conductivity for a light polarization par- 
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FIG. 2: Setups of an intelligent window, (a) VO2 on Si02. Reflectivity at high (dashed) and low temperatures (solid). Theory 
(red) : 60 nm VO2 on SiCh. Experiments : 50 nm VO2 on Si02 [261 ] (green), 50 nm VO2 on Pyrex glass [22] (blue), (b), 
(c) Setup with antireflexion coating : Ti02 / VO2 on SiC>2. Shown is the normalized visible transmittance (see text) an the 
corresponding color of the transmitted light for (b) high and (c) low temperature as a function of the layer thicknesses d V o 2 
and rf T io 2 - All theoretical data for E || [110] polarization, and a 3mm Si02 substrate. 



allel to the c-axis than for other directions [l2j. 

Having established the theoretical optical response of 
bulk VO2, and thus verified that our scheme can quan- 
titatively predict optical properties of correlated mate- 
rials, we now investigate the possibilities of VC>2-based 
intelligent window coatings [29|, [3(| . The effect to be ex- 
ploited can already be seen in the above bulk responses : 
The conductivities of both phases (Fig. Ufa), (b)) exhibit 
a close similarity in the range of visible light (w=1.7- 
3.0 eV), whereas in the infra-red regime (cj<1.7 eV) a 
pronounced switching occurs across the transition. As 
a result, heat radiation can pass at low external tem- 
peratures, while its transmission is hindered above T c . 
The insensitivity to temperature for visible light, in con- 
junction with the selectivity of the response to infrared 
radiation, is an essential feature of an intelligent window 
setup. Yet, for an applicable realization, other important 
requirements have to be met. First of all, the switch- 
ing of the window has to occur at a relevant, that is 
ambient, temperature. Also, the total transmittance of 
V02-films in the visible range needs improvement [29j ]. 
and the transmitted visible light should be uniform in 
order to provide a colorless vision. Experimentalists 
have addressed these issues and have proposed poten- 
tial solutions [H : Diverse dopings, M x V\- x O%, were 
proved to influence T c , with Tungsten (M=W) being 
the most efficient : A doping of only 6% results in T c w 
20°C [31j . However, this causes a deterioration of the 
infrared switching. Fluorine doping, on the other hand, 
improves on the switching properties, while also reduc- 
ing T c [32|, [33| . An increase in the overall visible trans- 
mittance can also be achieved without modifying the in- 
trinsic properties of the material itself, but by adding 
antireflexion coatings, for example TiC>2 [2tSj |. 

Here, we address the optical properties of window coat- 
ings from the theoretical perspective. In doing so, we 
assume that the specular response of VO2 layers is well- 
described by the optical properties of the bulk, and we 



use geometrical optics to deduce the properties of lay- 
ered structures. First, we consider the most simple of 
all setups, which consists of a single VC>2-layer (of thick- 
ness (2 v o 2 ) on a glass substrate [35j . Such a window has 
been experimentally investigated by Tazawa et al. j2?| 
and Jin et al. 26] . In Fig. [2^ we show their measured 
reflectivity data as a function of wavelength, in compar- 
ison with our theoretical results : At low temperatures 
(insulating VO2), the calculated reflectivity is in quan- 
titative agreement with the experimental data. In the 
visible range (A«400-700 nm), the reflectivity strongly 
depends on the wavelength. Therefore, the current win- 
dow will filter certain wavelengths more than others, re- 
sulting in an illumination of a certain color - an obvi- 
ous drawback. Moreover, the reflectivity in this region 
is rather elevated, causing poor global transmission. In 
the infrared regime (A>700 nm) - and beyond - the re- 
flectivity decreases, and radiation that causes greenhouse 
heating can pass the setup. At high temperatures, the 
infrared reflectivity switches to a rather elevated value, 
thus filtering heat radiation. The changes in the visi- 
ble region are less pronounced, but still perceptible, and 
both, the degree of transparency and the color change 
through the transition. The current setup is thus not yet 
suited for applications. 

We therefore now investigate a more complicated 
setup : An additional (rutile) Ti02-coating is added on 
top of the VQ2-layer, with the objective of serving as an 
antireflexion filter [26j. With the thicknesses g? V o 2 an d 
d Ti o 2 , the geometry of the current setup has two param- 
eters that can be used to optimize the desired optical 
properties. Since, however, each variation of them re- 
quires the production of a new individual sample under 
comparable deposition conditions, along with a careful 
structural characterization in order to guarantee that dif- 
ferences in the optical behavior are genuine and not re- 
lated to variations of the sample quality, the experimen- 
tal expenditure is tremendous. This led Jin et al. [26j to 
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first estimate a highly transmitting setup by using tabu- 
lated refractive indices and to produce and measure only 
one such sample. Here, we shall use our theoretical re- 
sults on VO2 to not only optimize the geometry (d TiQ2 , 
dvo 2 ) with respect to the total visible transmittance, but 
to also investigate the transmission color. Fig. [^displays 
the normalized visible specular transmittance [36| for the 
window in its high (b) and low (c) temperature state, as a 
function of both film thicknesses. On the same graph the 
resulting transmission colors are shown. The evolution 
of the light interferences in the layers causes pronounced 
changes in both, the overall transmittance and the color. 
The coating of VO2 globally degrades the transparency of 
the bare glass window. An increase of the Ti02-coating, 
on the other hand, has the potential to improve the total 
transmittance. This can be understood from the mech- 
anism of common quarter- wave filters. The wavelength- 
dependence of the real-part of the TiC>2 refractive in- 
dex, n Ti02 (A), results in an optimal quarter- wave thick- 
ness, <5 Ti02 (A)=A/(4ri TiO2 (A)), which varies from blue to 
red light only slightly from <5 TiO2 (A)=40 to 60 nm. This 
and the fact that the imaginary part of the refractive 
index, fc Ti02 (A), is negligible for visible light explains 
why the color does not change significantly with d Ti02 . 
While, as for TiC>2, the variation of the real-part of 
the VO2 refractive index yields a rather uniform ideal 
thickness <5 vo , (A), its imaginary part changes signifi- 
cantly (by a factor of 4) within the range of visible light. 
As a consequence, the color is very sensitive to VO2- 
deposition. At higher thickness d VQ2 , however, this de- 
pendence becomes smaller and the color lighter. Our 
theoretical transmittance profiles suggest relatively thick 
windows to yield good visual properties. Indeed, at low 
temperatures, Fig. [2jc), the local maximum that gives 
the thinnest window is located at (d Ti02 , e? V o 2 )~(40 nm, 
85 nm) within our calculation. However, this setup is 
still in the regime of important color oscillations. Given 
the uncertainties in industrial deposition techniques, it 
seems cumbersome to consistently stabilize colorless sam- 
ples. From this point of view, a thicker V02-film would 
be desirable. Indeed, while almost preserving the over- 
all transmittance, a colorless window at low temper- 
atures is realized for (d TiQ2 , <i VO2 )~(50 nm, 220 nm), 
or for (d T io 2 , d VO2 )—(>100 nm, 220 nm). In the high 
temperature state, Fig. [2fb), the transmittance is glob- 
ally lower than at low temperatures. Moreover, only 
(c? Ti o 2 , rfvo, )=(> 100 nm, 220 nm) yields a simultane- 
ously high transmittance in both states of the window. 

In conclusion, we have presented a novel scheme for 
the calculation of optical properties of correlated mate- 
rials, and applied it to vanadium dioxide, VO2. We find 
the bulk optical conductivity of both phases in quanti- 
tative agreement with experiments, and further validate 
our approach by comparing the transmittance of a VO2- 
layer on Si02 to experimental data. Finally, we optimize 
the geometry of a multi-layer setup of an intelligent win- 



dow, which uses the metal-insulator transition of VO2 to 
reduce the effect of greenhouse heating. This work can 
be considered as a proof of principle of the feasibility of 
theoretical materials design, since our techniques can be 
applied to the general class of correlated materials. 

We gratefully acknowledge discussions with T. Gacoin, 
G. Garry, A. Georges and L. Pourovskii. This work was 
supported by the Japan Science and Technology Agency 
(JST) under the CREST program, the French ANR un- 
der project CORRELMAT, and a computing grant by 
IDRIS Orsay under project number 081393. 



[2] 
[3] 
[1] 
[5] 

[6 

[r 
[8 
[9 
[10 

[11 
[12 

[13 

[11 
[15 
[16 
[17 



[18 

[19 
[20 
[21 

[22 

[23 
[24 

[25 

[26 
[27 

[28 
[29 

[30 
[31 

[32 
[33 



M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 

70, 1039 (1998). 

M. N. Baibich, et al, Phys. Rev. Lett. 61, 2472 (1988). 

G. Binasch, et al., Phys. Rev. B 39, 4828 (1989). 

Y. Paz, et al., J. Mater. Res. 10, 2842 (1995). 

C. G. Granqvist, Solar Energy Materials and Solar Cells 

91, 1529 (2007). 

G. Kotliar and D. Vollhardt, Physics Today 57, 53 
(2004). 

N. Bliimer, (PhD thesis, Universitat Augsburg, 2002). 
V. S. Oudovenko, et al., Phys. Rev. B 70, 125112 (2004). 
K. Haule, et al., Phys. Rev. Lett. 94, 036401 (2005). 
E. Pavarini, et al., New Journal of Physics 7, 188 (2005). 
L. Baldassarre, et al., Phys. Rev. B 77, 113107 (2008). 
J. M. Tomczak, (PhD thesis, Ecole Polytechnique, 2007). 

H. W. Verleur, A. S. Barker, and C. N. Berglund, Phys. 
Rev. 172, 788 (1968). 

K. Okazaki, et al., Phys. Rev. B 73, 165116 (2006). 
M. M. Qazilbash, et al., Phys. Rev. B 74, 205118 (2006). 
A. Khurana, Phys. Rev. Lett. 64, 1990 (1990). 
A. J. Millis, in Strong Interactions in Low Dimen- 
sions (Physics and Chemistry of Materials with Low- 
Dimensional Structures, 2004), vol. 25, p. 195ff. 
P. Baum, D.-S. Yang, and A. H. Zewail, Science 318, 788 
(2007). 

M. M. Qazilbash, et al., Science 318, 1750 (2007). 
V. Eyert, Ann. Phys. (Leipzig) 11, 650 (2002). 
M. S. Laad, L. Craco, and E. Miiller-Hartmann, Euro- 
phys. Lett. 69, 984 (2005). 

A. Liebsch, H. Ishida, and G. Bihlmayer, Phys. Rev. B 

71, 085109 (2005). 

S. Biermann, et al., Phys. Rev. Lett. 94, 026404 (2005). 
J. M. Tomczak and S. Biermann, J. Phys.: Cond. Matter 
19, 365206 (2007). 

J. M. Tomczak, F. Aryasetiawan, and S. Biermann 
(2007), submitted to Phys. Rev. Lett.. la7Xiv:0704.0902l 
P. Jin, et al., Appl. Phys. A 77, 455 (2003). 
M. Tazawa, P. Jin, and S. Tanemura, Appl. Opt. 37, 
1858 (1998). 

J. B. Goodenough, J. Solid State Chem. 3, 490 (1971). 
S. M. Babulanam, et al., Sol. Energy Mater. 5, 347 
(1987). 

C. G. Granqvist, Thin Solid Films 193-194, 730 (1990). 
M. A. Sobhan, et al., Sol. Energy Mater. Sol. Cells 44, 
451 (1996). 

W. Burkhardt, et al., Thin Solis Films 402, 226 (2002). 
W. Burkhardt, et al., Thin Solis Films 345, 229 (1999). 



5 



[34] We employ the CIE 1964 conventions, along with the 
daylight illuminant D65, see e.g. K. Nassau. The Physics 
and Chemistry of Color: The Fifteen Causes of Color. 
(Wiley Series in Pure & Applied Optics, 2001). 

[35] We suppose quartz glass, Si02- All auxiliary refractive 
indices are taken from Handbook of Optical Constants of 
Solids by Edward D. Palik, Academic Press, 1985. 

I 36 ] floon™ d\S{\)T(\)/ J™ °™ dXS(X), with the spectrum 



of the light source S(X), and the transmittance T = 1 — R, 
R being the specular reflexion given by Fresnel's formu- 
lae. We thus neglect absorption due to inhomogeneities 
that lead to diffuse reflexion. This is justified for our ap- 
plications to windows. Also, VO2 has a glossy appearance 
and hence a preponderant specular response. 



